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Abstract 

A discrete total variation calculus with variable time steps is presented in this letter. Using this 
discrete variation calculus, we generalize Lee's discrete mechanics and derive variational symplectic- 
" " | ~ energy-momentum integrators by Kane, Marsden and Ortiz. The relationship among discrete total 

variation, Lee's discrete mechanics and Kane-Marsden-Ortiz's integrators is explored. 

' Keywords. Total variation, Discrete mechanics, Symplectic-energy-momentum integrators 

^ ; 1 Introduction 

In 1980's, Lee proposed an energy-preserving discrete mechanics with variable time steps by taking 
(discrete) time as dynamical variable [g, ||, ||. On the other hand, motivated by the symplectic property 
of Lagrangian mechanics, a version of discrete Lagrangian mechanics has been devoloped and variational 
integrators that preserve discrete symplectic two form have been obtained Jll], [l2|, [l5|, [l^, |Tij ]. However, 
variational integrators obtained in this way have fixed time steps and consequently they are not energy- 
C , preserving in general. 

Obviously, the energy-preserving discrete mechanics and variational integrator are more preferable 
since solutions of the Euler-Lagrange equations of conservative continuous systems are not only symplectic 
but also energy-preserving. In order to attain this goal, some discrete mechanics with discrete energy 
conservation and symplectic variational integrators are needed to study. Recently, Kane, Marsden and 
Ortiz have employed appropriate time steps to conserve a defined discrete energy and developed what they 
called symplectic-energy-momentum preserving variational integrators in |Io| . Although their approach 
is more or less related to Lee's discrete mechanics, but the discrete energy preserving condition is not 
\ derived by the variational principle. 

• The purpose of this letter is to generalize or improve these approaches as well as to explore the 

relation among discrete total variation, Lee's discrete mechanics and Kane-Marsden-Ortiz's integrators. 
We will present a discrete total variation calculus with variable time steps and a discrete mechanics 
that is discretely symplectic, energy preserving and has the correct continuous limit. In fact, this discrete 
variation calculus and mechanics are a generalization of Lee's discrete mechanics in symplectic-preserving 
sense and can directly derive the variational symplectic-energy-momentum integrators by Kane, Marsden 
and Ortiz. 

This letter is organized as follows. In the next section, we remind total variation calculus for con- 
tinuous mechanics. In section 3, we present a discrete total variation calculus with variable time steps 
and a discrete mechanics, derive Kane-Marsden-Ortiz's integrators and explore the relation among our 
approach, Lee's discrete mechanics and Kane-Marsden-Ortiz's approach. We finish with some conclusions 
and remarks in Section 4. 

Before ending this section, we recall very briefly the ordinary variational principle in Lagrangian 
mechanics for later use. Suppose Q denotes the extended configuration space with coordinates (t, q l ) and 
Q^> the first prolongation of Q with coordinates (i, q\ q 1 ) [ fL3| . Here t denotes time and q l , i = 1, 2, ■ ■ • , n 
denote the position. Consider a Lagrangian L : — > R. The corresponding action functional is 



On 
O 
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defined by 



S(q i (t)) = / L{t,q\t),q\t))dt, (1.1) 

J a 

where q l (t) is a C 2 curve in Q. 

Hamilton's principle seeks a curve q l {t) denoted by C h a with endpoints a and 6, for which the action 
functional S is stationary under variations of q l (t) with fixed endpoints. Let 



V = f(t,q) 



d_ 

dq l 



(1.2) 



be a vertical vector field on Q. Here g = (g 1 , • • ■ , q n ). By a vertical vector field we mean a vector field 
on Q not involving terms of form £(i, g) Jj- Namely, time £ does not undergo variation. 

Let F e be the flow of V, i.e., a one-parameter group of transformations on Q: F £ (t, q r ) = (t, q l ). 

t = t, (1.3) 

q l = f/(e,t,q), (1.4) 

where 



g\e,t,q)=f(t,q) := 5q l (t). 



(1.5) 



e=0 



In other words, the deformation (1.31.4) transforms the curve q l (t) into a family of curves q l (e,t) in 
Q denoted by C e h a which are determined by 

t = t, (1.6) 
«* = <7 i (e,t,g(t)). (1.7) 

Thus, we obtain a (sufficiently small) set of curves C^ a around C\. Corresponding to this set of curves 
there are a set of Lagrangian and action functionals 



*«/'(./)) Sufu. in = f Lmt,t),j.f{e,i))dl 



(1.8) 



Now, we can calculate the variation of S at q(t) as follows 

d 



SS = 



de 



dL d dL 
dq l dt dq l 



(1.9) 



For the fixed endpoints, <f) l (a,q(a)) — <fi l (b,q(b)) — 0, the requirement of Hamilton's principle, SS = 0, 
yields the Euler-Lagrangian equation for q(t) 



dL d dL 
dq l dt dq % 



0. 



(1.10) 



If we drop the requirement of 4> l (a, q(a )) = ^(b, q(b)) = 0, we can naturally obtain the Lagrangian 
one form on from the second term in (1.9): 



h = T7T-dq , 



(1.11) 



where dq % are dual to ^j, dq l {-Jj~) = dj. Furthermore, it can be proved that the solution of ( l.lOj ) 
preserves the Lagrangian two form 



lu l := d0 L . 



(1.12) 
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On the other hand, introducing the Euler-Lagrange one form 



the nilpotency of d leads to 



dE(q\q l ) + -UJ L = 0. 
dt 



(1.13) 



(1.14) 



Namely, the necessary and sufficient condition for symplectic structure preserving is that the Euler- 
Lagrange one form is closed || [?], ||. 

2 Total variation for Lagrangian mechanics 

Consider a general vector field on Q 



(2.1) 



Here q = (q 1 , ■ ■ ■ , q"). Let F e be the flow of V. The variations of (t, q l ) £ Q are described in such a way 

= (*,$*), (2.2) 



where 



with 



t = f(e,t,q), «P = fl*(e,t,g), 



d_ 
~de 



f{e,t,q)=S(t,q) :=St, 



e=0 



d 

~dc 



g%t,q)=<j>%q) ■= V 



(2.3) 



(2.4) 



The deformations (2.3) transform a curve q l {t) in Q denoted by C a into a set of curves <f(e,£) in Q 



denoted by C e5 , determined by 

t = f(e,t,q(t)), j i = g i (e,t,q{t)). (2.5) 
Before calculating the total variation of S, we introduce the first order prolongation of V denoted as pr 1 ^ 



l?r 1 V = Z(t,q)-+<l> i (t,q)— i + a* (t , q, q) — % , 

where 

a%q i q) = D t <j>%q)-^D^{t i q), 
where D t denotes the total derivative with respect to t. For example 



(2.6) 



(2.7) 



at 



For prolongation of vector field and the related formulae, we refer the reader to |Tj| . 
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Now we calculate the total variation of S straightforwardly: 



de 

_ d_ 
~ de 

d 

~ de 



S{q\e,t)) 



L(t,q i (e,i),-q l (e,t))dt 
dt 

L{t, q l (e,t),^(e,t))^-dt 



e=0 J a 
b d 



dt 



dt 



(t = f(e,t,q(t))) 



de 



e=0 



L(t,(?(e,t),j-q l (e,t))dt+ / L(t,q\t),q l (t))D t {;dt 



dt dq 1 dq 
dL d ( dL 
~dt + ~dt 

dL 



dt+ I LDtidt 
dt 



dL d dL 
dq 1 dt dq 1 



dq % 



(2.8) 



Here we have made use of (UA), (2.6), (2.7) and 



d_ 
~de 



dt d d 
dt dt de 



t = D£. 



If £(a,q(a)) = £(b,q(b)) — and (j) l (a,q(a)) = (j) l (b, q(b)) = 0, the requirement of 5S = yields the 
equation from £, the variation along the base manifold, i.e. the time, 



dL 



dL 



m + dt{w" t - L )^°' 

and the Euler-Lagrange equation from cj> 1 , the variation along the fibre, i.e. the configuration space, 

dL d dL 



(2.9) 



dq 1 dt dq 1 



0. 



(2.10) 



Here £ and (j? are regarded as independent components of the total variation. However, there is another 
decomposition for the independent components, i.e. the vertical and horizontal variations, see Remark 2 
below. 

If L does not depend on t explicitly, i.e., L is conservative, ^ = 0, then ( |2.9| ) becomes the energy 
conservation law 



dt 



H = 0, 



H := 



dL 
dq 



■?-L 



We expand the left-hand side of (2.11) and obtain 



dL 



dt V d(f ^ / ' f)nl At f)M I Q 



dL d dL 



dq 1 dt dq 1 



(2.11) 



(2.12) 



Thus, for a conservative L, energy conservation is a consequence of Euler-Lagrange equation. This agrees 
with Noether theorem, which states that the characteristic of an infinitesimal symmetry of the action 
functional S is that of a conservation law for the Euler-Lagrange equation. For a conservative L, is an 
infinitesimal symmetry of the action functional S, and its characteristic is — q l . From Noether theorem, 
there exits a corresponding conservation law in the characteristic form 



dL d dL 
dq 1 dt dq 1 



0. 



(2.13) 



4 



If we drop the requirement 



e(a, g(a)) = f (6, g(6)) - 0, ^(a, g(o)) = ^(b, q(b)) = 0, 
we can define the extended Lagrangian one form on from the second term in ( |2.S| ) 



(2.14) 



(2.15) 



Suppose g l (t,v q i) is a sol ution of ( 2.10 ) depending on the initial condition v q i £ Restricting q l (e,t) 

to the solution space of ( 2.1C ) and using the same method in pijl , it can be proved that the extended 
symplectic two form is preserved 

(prVTft L = Q L , Q L := d§ L . (2.16) 
where pr 1 g l (s, v q i)) = (s, g l (s, v q i), j^g l (s,v q i)) denotes the first order prolongation of g l {s,v q i) |T^| . 



Remark 1. If £ in (2.1) is independent of q, the deformations in (2.3) are called fiber-preserving. In this 
case, the domain of definition of q l (e,t) only depends on the deformations (2.3). While in the general 
case, the domain of definition of tf*(e, i) not only depends on the deformations (2.3) but also on q r {t). 



Remark 2. Using the identity 

dL 
dt 



dL .. 
37 1 tttV 



dt \dq 



L 



dL 
dq l 



d dL 

didq 7 1 9 '' 



the Eq. ( p.8|) becomes 
SS = 



(2.17) 



(2.18) 



According to (2.4), <p l — Sq l should be regarded as the total variation of q\ 5q r = 8yq l + Suq 1 , since 
the variation of t also induces the variation of q l denoted as Snq 1 , the horizontal variation of q l . Due to 
£ = St in (2.4), the horizontal variation of q % should be 8uq % — £q l and consequently 4> l — £q l is interpreted 
as vertical variation 5yq\ i.e. the variation of q l (t) at the moment t (see, for example |jj). Therefore, 
the first two terms in ( 2.18| ) come from vertical variation Syq 1 and the last term comes from horizontal 
variation St. The h orizon tal variation of S with respect to the horizontal variation Snq 1 = S.q 1 just gives 
rise to the identity ( 2.17 ). 



3 Discrete mechanics and variational integrators with variable 
time steps 

In this section, by means of a calculus of discrete total variations we develop a discrete Lagrangian 
mechanics, which includes the boundary terms in Lee's discrete mechanics that give rise to the discrete 
version of symplectic preserving. The discrete variation calculus is mainly analog to Lee's idea that 
(discrete) time is regarded as dynamical variable, i.e. the time steps are variable \ % p[ p| . And the 
vertical part of this discrete variation calculus is similar to the one in [|l0| [ll], |l2|, |l5|, lq, |17|. Using this 
calculus for discrete total variations we naturally derive Kane-Marsden-Ortiz's integrators. 

We use Q x Q for the discrete version of the first prolongation for the extended configuration space Q. 
A point (to, ffo; ii) qi) £ Q X Q Corresponds to a tangent vector . A discrete Lagrangian is defined 
to be L : Q x Q — > K and the corresponding action to be 



N-l 



, <Zfc, tk+i, <7fc+i)(ifc+i — tfc). 



(3.1) 



k=0 



1 In this section q is an abbreviation of (q 1 , q 2 , ■ ■ ■ ,q n ). 
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The discrete variational principle in total variation is to extremize § for variations of both q k and t k 
with fixed endpoints (to,9o) and (tjVi Qn)- This discrete variational principle determines a discrete flow 
$:QxQ->QxQby 



$(4-1) <7fc-i>4, qk) — (4, 9fc, 4+i, ?a+i)- 



(3.2) 



Here (4+i ; <?&+i) are found from the following discrete Euler-Lagrange equation, i.e. the variational 
integrator, and the discrete energy conservation law (for conservative L) 



and 



(4+1 ~ tk)D2h(tk, Qk>tk+1, Qk+l) + (4 ~ tk-l)D4h(t k _i,q k _i,t k ,q k ) — 0, 

{t k+ i - t k )D{L(t k ,q k ,t k+ i,q k+1 ) + D s la(t k -t, q k -i, t k , q k ){t k - 4-i) 
-h(t k ,q k ,t k+ i,q k+ i) + L(t k -i,q k -i,t k ,q k ) = 0, 



(3.3) 



(3.4) 



for all k £ {1, 2, ■ • • , N— 1}. Here Di denotes the partial derivative of L with respect to the ith argument. 
The Eq. (3.3) is the discrete Euler-Lagrange equation. The Eq. (3.4) is the discrete energy conservation 
law for a conservative L. The integrator ( 3.3 )-( 3.4 ) is just Kane-Marsden-Ortiz's integrator. 
Using the discrete flow the Eqs. (3.3) and (3.4) become respectively 



{t k+ i - t k )D 2 h o $ + (t k - t k ^)Dih = 0, 
((tk+i - tkjDtL - L) o $ + D 3 L + L = 0. 



(3.5) 
(3.6) 



If (4+1 — t k )D2L and (t/c+i — tk)D\L — L are invertible, the Eqs. ( |3.5| ) and (3.6) determine the discrete 
flow $ under the consistency condition 



((4+1 - * fc )£>iL - L)- 1 o (£> 3 L +L) = - ifc)^)" 1 o (4 - 4-1)1)41 



(3.7) 



Now we prove that the discrete flow <£> preserves a discrete version of the extended Lagrange two form 
Ql- As in continuous case, we calculate dS for variations with varied endpoints. 

d§(io,go, • • • ,t N ,q N ) ■ (St ,Sq , ■ ■ ■ ,5t N ,6q N ) 

JV-l 



fc=0 
JV-l 



E {D 2 L(t k ,q k ,t k+1 ,q k+1 )Sq k + D 4 L(t k ,q k ,t k+1 ,q k+1 )Sq k+1 )(t k+1 - t k ) 

k=0 
N-l 

, gfc, t k +i,q k +i)5t k + D?,L{t kl q k ,t k+ i,q k+ i)St k+ i)(t k+ i — t k ) 

fc=0 
(V-1 

E L(t kl q k ,t k+1 ,q k+1 )(6t k+ i — 8t k ) 
k=o 
-1 

E ^(4 ; %,4+i) <?fc+i)(4+i — tk)Sq k 
=0 

N 

y^DjL^k-i, g k -i, t k , q k ){t k - t k -i)Sq k 
fe=i 

iV-l 

, ?fc)4+i)3fe+i)(4+i — t k )St k 

fc=0 
AT 

E^3i(4-i , %-i,4,%)(4 - tk-i)St k 
fe=i 

TV — 1 AT 

E ^(4, ?fe, 4+i, ?fc+i)(-^4) + E £(4-i, gfe-i) 4, gfc)^4 



fe=0 
AT 



fe=0 



fc=l 
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N-l 

- J~] (D 2 L(t k , gfc, tfc+i, qk+i)(tk+i - tk) + DiL(t k -i,qk-i,tk,qk){tk - *fe-i))^gfc 

k=l 
N-l 

+ ^ {D\L(t kl q k , tk+i,qk+i){tk+i — tk) + £>3£(ifc-i, <Zfc-i, ifc, qk){tk ~ t k -i) 

fe_1 + L(t k _ 1 ,q k _ 1 ,t k ,q k ) — L(t k , q k ,tk+i, <7fc+i))£ifc 
+ L> 2 i(^o,go,ii,9i)(*i ~t )5q + -D 4 -^(*JV-i, &v)(*iv - t N -!)8q N 
+ (L>iL(i ,<7o,ii,gi)(ii - £o) - L(*o, go, 9i))^o 

+ {D3L(tN-l,qN-l,tN,qN){tN — tN-l) + L(tff-i, QjV-lj *iV, <Ziv))<ftiV- 



(3.8) 



We can see that the last four terms in ([3^) come from the boundary variations. Based on the boundary 
variations, we can define two one forms on Q x Q 



and 



= (Dih(tk,qk,tk+l,qk+l)(tk+l — tk) — ^(tk,qk,tk+l,qk+l))dtk 

+ D 2 h(tk,qk,tk+i,qk+i)(tk+i — t k )dq k , 



\ (tk,Qk,tk+l,qk+l) 

— (D^Litk, q k ,t k+ i, q k+ i)(t k+ i — t k ) + l^(t k ,q k ,t k+ i,q k+ i))dt k+ i 
+ D4h{tk, q k , t k+ i,q k +i)(t k+ i — t k )dq k+ i. 



(3.9) 



(3.10) 



Here we have employed the notations in pi)]. W e regard the pair , 9^) as the discrete version of the 
extended Lagrange one form i9l defined in Q2.15 ). 

Now we parameterize the solutions of the discrete variational principle by the initial condition (to, qo,ti,qi) 
and restrict § to that solution space. Then Eq. (|3.8|) becomes 



(3-11) 



d^(t ,qo, ■■ ■ , tN,qN) ■ (St , Sq , • • • , St N , Sq N ) 
= Qi{to,qo,h,qi) ■ (^0-^0,^1,^1) 

+ #l {tN-l, qN-1, tN, 97v) ■ (8tN-l,5qN, 5tN-i,5qN-i) 

= 0£(*o,?o,*i,9i) • (Sto, 6q , Stx,8qi) 

+ ($ w - 1 )^+(t ,g ,ii,9i) • (^0,^0,^1,^1). 

From ( fo.ll| ), we can obtain 



The Eq. (3.12) holds for arbitrary N > 1. By taking N=2, it leads to 



ds = e- + $*e+. 



Taking exterior differentiation of (3.13), it follows that 

<r«) = -de h 



From the definition of L and 9t, we know that 



d{Ut k +i-t k )). 



Taking exterior differentiation of (3.15), we obtain d9^ = —d9 h ■ Define 

S^L = d6^ = — d9, . 



(3.12) 



(3.13) 



(3.14) 



(3.15) 



(3.16) 
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Finally, we have shown that the discrete flow $ preserves the discrete extended Lagrange two form fi^. 

$*(fi L ) = n L . (3.17) 



Now we show that the variational integrator (3.3), the discrete energy conservation law (3.4) and the 
discrete extended Lagrange two form Ol converge to their continuous counterparts as tk+i — * tk, t k -i — * 



Consider a conservative Lagrangian L(q, q). For simplicity, we choose the discrete Lagrangian as 



L(tfe, qk,tk+i, 3fe+i) 
The variational integrator (ET3) becomes 



gfc+i — gfc 

tk+l ~ tk 



dL 
dq k 



(q kl A t q k ) 



1 



where A t q k 



q k+ l~q k 



Atqk 



dL 

tk+i — tk \dA t q k 

qk—Qk-i 



(q k ,A t q k ) 



dL 



dA t q k - 



■(q k -i,A t qk-i) = 0, 



(3.18) 



(3.19) 



It is easy to see that as tk+i 



tk, tk-i 



tk, the Eq. ( pMS| ) converges to 
0. 



dL d dL 
dq k dt dq k 



The discrete energy conservation law (| 



J) becomes 

Ek+i — Ek 
tk+i — t k 



= 0, 



(3.20) 



(3.21) 



where 



E k +i 
E k = 



-A t q k - L[q k , 



dA t q k £fc+i — ifc 

dL qk-Qk-i, 
-A t q k -i - L{qk-i, )• 



dA t q k -i 



tk — t k -l 



The Eq. (3.21) converges to 



d f dL 



dt \dq k 



qk 



L 



(3.22) 



as tk+i — » tk, tk-i tk- 

N ow w e consider the discrete extended Lagrange two form de fined by ( 3.16 ). Under the discretiza- 
tion (|3.18 ), the discrete extended Lagrange one form 9^ defined in ( 3.10) ) becomes 



9+ 



L(q k ,A t q k ) 



dL 
dA t q k 



A t qk dt k 



dL 
dAtqk 



dqk 



+ 1- 



(3.23) 



From ( 3.23j ), we can deduce that 9^ converges to the continuous Lagrangian one form defined by 
( 2.15 ) as tk+i — > tk, tk-i — ■> tk- Thus, we obtain 

fi L = d9+ -» d& L = fl L , (t k+1 -» t h , tk-i -» t fc ). (3.24) 

The variational integrator ( ^.3[ ) with fixed time steps does not conserve the discrete energy exactly in 
general but the computed energy will not have secular variation [||, Q . In some cases such as in discrete 
mechanics proposed by Lee in [||, [| [I] , the integrator ( [3~3] ) is required to conserve the discrete energy 
(3.4) by varying the time steps. In other words, the time steps can be chosen according to (3.4) so as 
to the integrator (3.3) conserves the discrete energy (3.4). The resulting integrator also conserves the 
discrete extended Lagrange two form d6£ . This fact had not been discussed in Lee's discrete mechanics. 
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4 Conclusions 



We have presented the calculus of total variation problem for discrete mechanics with variable time 
steps referring the one for continuous mechanics in this letter. Using the calculus for discrete total 
variations, we have proved that Lee's discrete mechanics is symplectic and derived Kane-Marsden-Ortiz's 
integrators. It is well acknowledged, an energy-preserving variational integrator is a more preferable and 
natural candidate of approximations for conservative Euler-Lagrangian equation since the solution of 
conservative Euler-Lagrangian equation is not only symplectic but also energy-preserving. 

As is mentioned, Kane-Marsden-Ortiz's integrators are closely related to the discrete mechanics pro- 
posed by Lee ||, f|. In Lee's discrete mechanics, the difference equations are just Kane-Marsden- 
Ortiz's integrators. However, Lee's difference equations are solved as boundary value problems, while 
Kane-Marsden-Ortiz's integrators are solved as initial value problems. 

Finally, it should be mentioned that two of the authors (HYG and KW) and their collaborators have 
presented a difference discrete variational calculus and the discrete version of Euler-Lagrange cohomology 
for vertical variation problems in both Lagrangian and Hamiltonian formalism for discrete mechanics 
and field theory very recently g|. In their approach, the difference operator with fixed step- length 
is regarded as an entire geometric object. The advantages of this approach have already be seen in the 



last section in the course of taking continuous limits although the difference operator A t in ( 3.1q ) is 
of variable step-length. This approach should be able to be generalized for the discrete total variation 
problems. 
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